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A direct method to reconstruct the transverse structure of the baryon source 

formed in a relativistic heavy ion collision is presented. The procedure makes 

CN . use of experimentally measured proton and deuteron spectra and assumes that 

, deuterons are formed via two-nucleon coalescence. The transverse density shape 

, and flow profile are reconstructed for Pb+Pb collisions at the CERN-SPS. The 

O . 

Q\ . ambiguity with respect to the source temperature is demonstrated and possible 

ways to resolve it are discussed. 
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■ One of the main goals of modern nuclear physics is to understand the behavior 
of nuclear matter under extreme conditions. This can provide a solid basis for a 
reliable description of the nuclear equation of state and therefore answer many unsolved 
problems in astrophysics and early universe cosmology. In laboratory experiments there 
is only one way towards this goal, which is to study high energy collisions between heavy 
nuclei. 

The complicated non-equilibrium dynamics associated with such collisions are nowa- 
days the object of intensive theoretical investigations. There have been many attempts 
to develop a realistic approach describing the full space-time evolution of the system 
0, 0. Furthermore, due to the large number of produced particles and their sub- 
sequent re-scattering, a certain degree of local thermodynamic equilibrium might be 
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established at an intermediate stage of the reaction [|, ^ ||. One can therefore try to 
use simpler models, dealing with the properly parametrised phase space distribution 
functions for different particle species. 

However, on this road one always encounters an inherent problem. What can be 
measured in an experiment are only the momentum spectra of particles. Each measure- 
ment is only a projection of the full phase space information on the space-time points 
where particles decouple, the so-called freeze-out region. The fundamental ambiguity 
that several phase space distributions, after summing over the freeze-out coordinates of 
particles, can lead to the same single particle spectrum is left, in principle, unresolved. 
We see that one must disentangle what can be referred to as an inversion problem. 

Although it is clear that independent single particle spectra are not sufficient to 
resolve the ambiguity, we fortunately have at our disposal a promising probe, the 
deuteron (light nuclear clusters, in general). Previous studies have shown that 
deuteron production cross sections can be understood in terms of the phenomenological 
coalescence model. It is founded on the fact that a neutron and a proton can fuse into 
a deuteron only if their relative distance in position and momentum space is small, on 
the order of the size of the deuteron wave function. Furthermore, because of the small 
binding energy (~ 2 MeV), deuterons can survive break-up only when scatterings are 
rare. These two constraints reveal that deuteron production can only take place at 
freeze-out. This is the reason why, with a suitable use of single and composite particle 
spectra, one can constrain the freeze-out phase space distribution of protons. 



Another important observation made in recent experiments @, rU| is the devel- 
opment of collective effects during the reaction. The collective behavior is revealed by 
increasing inverse slopes of spectra for particles of increasing mass. This can be un- 
derstood in terms of collective flow of nuclear matter, leading to a position-dependent 
local velocity. 



In this letter we extend our previous work [0], attempting to reconstruct the phase 
space distribution of protons directly from the observed proton and deuteron spectra, 
exploiting the coalescence prescription together with the notion of collective flow. It 
is assumed that one can characterise the system at freeze out by a position-dependent 
collective velocity and a position-independent local temperature Tq. We are mainly 
interested in the transverse dynamics at mid-rapidity for central collisions of large 
symmetric systems. We therefore make use of a relativistic description of collective 
flow based on the Bjorken picture for the longitudinal expansion, together with a 
longitudinally-independent transverse velocity v±{r±). The contraction of the four- 
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momentum with the collective four-velocity can be written as |12 



Pf,u''{x) =-f±{r_i) (m_L cosh(?/ - i]) - p_i_- t;j_(rj_)) , (1) 



where m_|_ = y + is the transverse particle mass while y and r] are the momentum 
and space-time rapidities, respectively. Choosing a transverse profile np{r±) for the 
local density, we can write the proton phase space distribution as 

fp{x,p) = (27r)3e-^A^"^^^)/^0i?pnp(rx). (2) 

From now on we drop the subscripts in 7^, v± and r±, understanding that these quanti- 
ties denote the transverse degrees of freedom. In the above expression the normalisation 
constant for the Boltzmann factor in the local frame is defined as 

B;' =jd'p e- y = 4vr m^T, K,{m/T,) , (3) 

where K2 is the modified Bessel function of second order. The local density must 
be normalised to the measured differential multiplicity dNp/dy. To do this we first 
calculate the invariant momentum spectrum. It can be obtained by integrating the 
phase space distribution on the freeze-out hypersurface, using the Cooper-Frye formula 

d^N 1 r 



dyd^p^_ (27r)3, 

The form of the integration measure follows from our simplifying choice of freeze- 
out hypersurface as a sheet of constant longitudinal proper time r = tq. We neglect 
longitudinal edge effects which are negligible for the mid-rapidity region of the spectra. 
We also assume a simultaneous freeze-out in the transverse direction and neglect surface 
emission at r < tq. Within this assumption the integration measure assumes the form 

p'^da^ = Tq m±_ coshiy — rj) d^rdrj . (5) 

After substituting this expression in eq. (^, the 77 integration can be done analytically. 
Now the rapidity density can be obtained by integrating out the transverse momentum 
dependence of Sp{pi_). Finally, one arrives at the normalisation condition for the local 
proton density 

—-^ = 2tctq drr 7(r) np(r) . (6) 
dy J 

This completes the definition of the proton phase space distribution. 

The deuteron phase space distribution can be calculated on the basis of the coales- 
cence model. In the density matrix formalism [13] one obtains 

fdix,p) = lj ^j^^ fp{x+,p+) f nix-, p^) Pdiy,q} , (7) 
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where the factor | accounts for the spin-isospin couphng of the neutron-proton pair 
into a deuteron state. It is assumed that neutrons (not measured) evolve in the same 
way as protons and that fn = Rnp fp, where i?„p = 1 — 1.5 is the neutron to proton ratio 
in the source. In the following we will take Rnp = 1.2. The phase space coordinates of 
the coalescing pair are x± = x±y/2 and p± = p/2±q while is the Wigner density for 
the deuteron relative motion. The coalescence prescription is greatly simplified when 
considering large and hot systems. In this case one can neglect the smearing effect of 
the deuteron Wigner density in comparison to the characteristic scales of the system 
in position and momentum space |14|. Then the deuteron phase space distribution 



becomes 

fd{x,p)^lRnp[fp{x,p/2)f . (8) 

This expression simply means that two nucleons, each with 4-momentum p/2, form a 
deuteron with 4-momentum p at space-time point x. 

Repeating the same reasoning for deuterons as done above for protons and inserting 
the corresponding expressions for fp and fd into eq. (^), one obtains the relation 

nd{r) = \dnl{r) (9) 

between the local densities of deuterons and protons. The proportionality coefficient 
has dimension Lr" and carries information on the characteristic scales in the problem. 
Its explicit form is 



A, = :-Rnp{2ixf^ . (10) 

We now have built up the framework to establish the reconstruction procedure. In 
a fashion which holds both for protons and for deuterons, using eqs. (j^), and (j^), we 
can calculate their invariant momentum spectra. Due to our simple choice of freeze-out 
hypersurface the integrations over the space-time rapidity and the azimuthal angle can 
be easily performed, leading to the expression Q 

S{p±)=C drrKi{ — ) Jo( )ron(r), (11) 

J Jo Jo 

where C = ATTBm± and Ki and Jq are modified Bessel functions of first and zeroth 
order. Here one can explicitly see the ambiguity in the description of the individual 
single particle spectrum. The two functions v{r) and n{r) cannot be mapped out 
uniquely from only one function as with the transverse momentum spectrum. This is 
true if protons and deuterons are treated independently, but with the link provided 
by the coalescence model, this ambiguity can, at least partially, be removed. Let us 
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introduce for convenience a new, dimensionless, density function 

n{v) = ^-^Ton{r) (12) 
V av 

and change the integration variable in eq. (0) from r to f . In this way we obtain the 
new expression for the momentum spectrum 

S{p^) = c[dvv K,C^) /o(^) fi{v) , (13) 

from which the one-to-one correspondence between h{v) and S{p±) is evident, since all 
the functions in the integrand are single-valued. The normalisation condition (^) for 
np(r), due to the definition (p!2D, now becomes 

dN r 

—^ = 2TTdvv-fnJv), (14) 
dy J 

together with the analogous one for n^. It is instructive to consider the limit Tq — in 
eq. (p!3|). Using the asymptotic expression for large argument for the Bessel functions Ki 
and Jo and performing the integration with the saddle point method, one obtains that 
V = p±/m± and m± = m'j. The transverse momentum spectrum is simply expressed 
through h as 

TTl 

S{p±) ^n{p±/mj_) . (15) 

One can easily find n when the observed spectrum has an exponential shape with 
inverse slope as 

(mT*)'^/^ 

Expressing the momentum variable in terms of velocity, p± = mvj, we obtain 

niv)-l'bfe-^^^^-^\ (17) 

with bo = m/T^. Observe that the exponent behaves as exp (— 6o "^^^72) for small v 
and as exp (— bo 7) for large v. In the limiting case Tq ^ the function n is therefore 
uniquely determined from the spectrum. A finite temperature introduces an inevitable 
intrinsic smearing, which precludes an exact determination of n for a single particle 
species. 

The role of the coalescence model is now to establish the link between proton and 
deuteron spectra. Substituting in eq. (H) the definition of h from eq. ([T^) , both for 
protons and for deuterons, we obtain the simple differential equation 



n,{v) = - ^ n^{v) , (18) 
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which can be directly integrated and leads to the closed solution 



r' = 2— / duu^^. (19) 
To Jo nd{u) 

Therefore, by independently extracting the functions Up and from the observed 
momentum spectra Sp{p±) and Sd{p±), we can find the function r{v) by a simple 
numerical integration. Inverting the obtained function as r{v) v{r), we obtain the 



collective velocity profile. Substituting it in eq. (|T3), we finally obtain the local proton 
density 

^vjrldv^ npivjr)) 
r dr tq 

We now explore more closely the content of eqs. (|T^ and (^). Taking the limit for 
small transverse velocity in eq. (|T^, one finds that v{r) c:^ H r and the dimensionful 
scale in position space is set by the "Hubble constant" 

The r.h.s. of this expression is obtained from eqs. @ and (0), using the asymptotic 
expression for large argument for the Bessel function K2 in the Boltzmann factors. 
Using this result and taking the limit for r ^ in eq. (pOD , we obtain the proton 
spatial density at the origin 

SKnp V TT / np[{}j 

This result will be used in the following discussion. 

The described procedure to reconstruct the proton phase space distribution can now 
be applied to the transverse momentum spectra, measured by the NA44 collaboration 



at the CERN-SPS, for Pb+Pb collisions at 158 A GeV [|T5[. Using these data, we fitted 
rip and ha, using eq. (|^) for protons and deuterons and assuming different values for Tq. 
Guided by eq. (|1^) we chose a rather flexible form of the profile function, characterized 
by three parameters, k, a and b, as 

n[v) ~ 7 e ^ ' ^ ^ ' ' . (23) 

The exponential factor has now the limiting behavior exp(— af^) for small f, while 
exp i^—b'y) for large v. We separated the scales a and b for small and large f , in order 
to have more freedom in fitting the curvature of the spectrum. The parameters fcp. 
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Op, bp, kdi ad, hd, are extracted from the experimental data with a Monte Carlo search 
minimising 

The only constraint imposed was 6d > 26p: because of our choice of h in eq. (|23|) , the 
integrand in eq. ([19D is exponentially divergent for v 1, thereby giving the limit 
v{r) — > 1 for r oo. The values obtained in this way are listed in Tab.^ and the 
reconstructed function rip is shown in Fig. |l]. We do not show the deuteron profile since 
it is very similar to the proton one. The fitted spectra, normalised as dNp/dy = 22 
and dNd/dy = 0.3, are shown in Fig. ^ One clearly sees that the calculated spectra for 
different temperatures are indistinguishable from one another. On the other hand, the 
reconstructed profiles n are very different for different temperatures. They are very 
much peaked at large temperatures and become broad as the temperature drops. This 
is simple to understand, since a lower temperature introduces less momentum spread 
than a higher one and n compensates for that by broadening. One can notice that 



the high-temperature profiles have a structure resembling a blast wave ||T6[ where most 
particles have approximately the same velocity. 

In all calculations the freeze-out time was fixed at tq = 10 fm/c After numer- 
ical integration of eq. (|19|) , we obtained the function v{r) shown in Fig.^ for different 
temperatures. It shows a linear rise at small r and saturates for large r. The velocity 
profiles clearly depend on the temperature chosen. 

Making use of eq. (|20|) we also obtained the local proton density, plotted in Fig.§. 
One can observe a behavior similar to fip. At high temperature the density shows a 
shell-like structure which disappears as the temperature is lower. It should be em- 
phasised that the shapes of v{r) and np{r) are quite sensitive to the curvature of the 
momentum spectra. 

From the plots of v{r) in Fig. ^, one might have the impression that flow is somehow 
stronger at high temperature. This is not the case, since the apparent flattening of 
v{r) at small r is accompanied by a higher saturation value and a corresponding larger 
extension of np{r). To see this more quantitatively we calculated the mean-square 
value for the velocity field. The mean, taken with respect to the proton density in the 
global frame defined by pp{r) = 7(r) np{r), is calculated as 

(^')= f^l '[d'rv\r)ToPp{r), (25) 



dy 

and the results are listed in Tab. |l] for different temperatures. As one can see, (f^)^/^ 
indeed grows with decreasing Tq. 
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The broadening of Hp as the temperature decreases is dictated by the conservation 
of the total phase space volume. The particular value chosen for tq also affects the 
resulting transverse extension of the source. All of this is evident after examining 
eq. (|2T|). A large temperature results in a large H, i.e. in a smaller transverse extension 
scale H^^. As a consequence the average proton density is larger. This is explicit in 
eq. (P2|). To quantify the transverse size of the source we have also calculated the 
mean-square transverse radius 

{r') = (^] 'fd'rr'roppir). (26) 



dy 

The resulting values of (r^)^/^ are listed in Tab.jl]. One can observe the increasing 
mean transverse size with decreasing temperature. It is therefore necessary to know Tq 
precisely in order to determine the system size. This information cannot be extracted 
solely from proton and deuteron spectra. Contrary to what is commonly done, source 
radii cannot be extracted from the d/p^ ratio (usually assumed to be inversely pro- 
portional to the source volume) unless Tq is known or further assumptions are made. 
One possibility relies on a rough estimate of the characteristic value of Up at freeze-out. 
Taking the total hadron number density to be ~ Po/S, where po = 0.17 fm~^ is the 
nuclear saturation density, with an average hadron-hadron cross section of ~ 30 mb, 
we have a mean free path of ~ 6 fm. This is long enough to reach the freeze-out 
conditions, especially for a rapidly expanding system. Since the ratio of all hadrons 
[p, n, n, K, . . .) to protons is ~ 15 for Pb+Pb collisions at 158 AGeV |jl8|, we can 
estimate that 

11 ^ Q 

n„ ~ — X - X 0.17 fm~^ ^ 0.0038 fm"^ . (27) 
15 3 

With this number, one can say that our reconstruction procedure suggests the preferred 
values of Tq ~ 100 MeV and {r^y/^ - 9.2 fm, as follows from Fig.| and Tab.|l|. 
This is also consistent with the temperature dependence in eq. (^2]). Comparison with 
microscopic simulations of the freeze-out distributions at AGS energies [|19|, ^ and at 
SPS energies PT| also shows that the temperature around 100 MeV would be preferable. 
To make a definite statement, one must finally keep in mind that the average source 
size is influenced by the choice of tq. A larger (smaller) value would result in a smaller 
(larger) average size. Internal consistency requires that tq is large enough so that the 
flow fleld has time to develop and the source to expand transversely during the reaction. 
Details about these isuues can only be discussed within a dynamical calculation, unless 
an independent measurement of tq is available. 

It is clear that to resolve the remaining ambiguity in the source temperature one 
needs some additional experimental information. Recently, tt tt HBT correlations data 
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were used as a constraint |2^. Since pions may freeze-out in a different way than 
protons, it would be even better to consider pp HBT correlations, although they are 
much more sensitive to final state interactions than pions. We would also like to 
mention that the neutron phase space distribution, here assumed to be proportional to 
the proton one, could be also reconstructed with a similar procedure for triton spectra. 
Furthermore, heavier clusters are more sensitive to collective flow than to temperature 
2m , providing additional constraints. 



In conclusion, we have demonstrated how the proton phase space distribution at 
freeze-out can be reconstructed from the measured transverse momentum spectra of 
protons and deuterons. The calculations, made with several simplifying assumptions, 
show that the proposed method gives meaningful results, although the ambiguity with 
respect to the temperature of the source remains. Two modifications will make this 
approach more realistic. First, the deuteron size should be explicitly included in the 
calculation. The wave function's tail plays a quantitatively important role, especially 
at high energies, when the large number of produced particles forces nucleons to be 
far apart at freeze-out. Second, the freeze-out hypersurface itself is more complicated 
than a constant proper time hyperbola. A correct description of early evaporation of 
fast particles is needed in order to have the high momentum part of the spectra under 
control. Nevertheless the described procedure provides a basis for more quantitative 
analyses, possibly including HBT correlations for protons. 
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tive suggestions and Ian Bearden, Jens J0rgen Gaard0je, Allan Hansen and the NA44 
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Figure 1: Profile rzp, characterized by the parameters shown in tab.|l[ 
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Figure 2: Transverse momentum spectra for protons (top) and deuterons (bottom), 
plotted as functions of m_|_ — m. They are obtained by fitting the NA44 data ||15|| , 
using the trial functions of eq. (|T3D with the parameters shown in Tab. |l[ Calculated 
curves corresponding to different temperatures are superimposed and indistinguishable. 
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Figure 3: Reconstructed transverse flow velocity flelds for different temperatures. The 
proflles shows a strong rise at small r, leveUing off at large r. 
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Figure 4: Reconstructed local density of protons as a function of the transverse coor- 
dinate. The resulting profile shows a shell-like structure at high temperature, turning 
into a more uniform behavior at smaller Tq. 
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Table 1: Fitted parameters characterizing the profiles hp and hd for different values of 
the source temperature Tq. Root-mean-square values for velocity and transverse radius 
are also hsted. 



To 
(MeV) 


kp 


ttp 


bp 


(^2^1/2 


(^2^1/2 

(fm) 


kd 


ad 


bd 


50 


91.56 


47.63 


52.40 


0.54 


17.60 


165.29 


83.55 


107.77 


75 


108.50 


51.37 


68.91 


0.51 


11.93 


214.21 


106.65 


146.04 


100 


75.86 


34.19 


66.45 


0.48 


9.20 


146.62 


67.85 


133.62 


125 


68.48 


23.18 


105.00 


0.44 


8.49 


67.24 


9.78 


211.05 


150 


172.89 


17.47 


602.83 


0.39 


7.02 


307.96 


31.06 


1211.85 
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